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Abstract 
Compared to the gigahertz (GHz) permeability of single Fe nanowire, terahertz (THz) 


permeability of single L1o-FePt nanowire at its remanent state (Mr/Ms =1.0) has been simulated in 
the micromagnetics approach. Due to its very large magnetocrystalline anisotropic field, it is 
feasible to obtain strong THz permeability at 0.348 THz , which are ever thought only possible in 
some metamaterials. The THz permeability spectra of Llo-FePt are shown obviously different 
from those of Fe nanowire at gigahertz (GHz). Unusual negative imaginary parts of permeability 
(u"«0) is found related to the equivalent negative damping constant, which is explained from the 


perspective of abnormal precession of natural resonance at THz for FePt nanowire . 
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1. Introduction 


The magnetic permeability characterizes the response of magnetic moments in magnetic 
materials to external magnetic fields. For electromagnetic applications, it is crucial to understand 
how the complex permeability values (real parts and imaginary parts) vary with the frequency [1]. 
The resonance frequency (f;) of permeability spectrum governs the operation scope of magnetic 
materials. For applications demanding low magnetic loss, the resonance frequency should be far 
away from the working frequency, such as magnetic cores of transformers, microwave magnetic 
devices (such as circulators). For applications demanding large magnetic loss, the resonance 
frequency should fall into the operation scope so as to dissipate the electromagnetic wave energy 
via the magnetic resonance phenomena, for instance, the anti electromagnetic interference or 
stealth technology. Up to now, it is believed that the highest operating frequencies of 
electromagnetic devices employed naturally synthesized magnetic materials falls into the W-band 
(75-100 GHz), which was realized by employing BaM type (BaFe:2019) hexaferrites in which Fe 
cations are properly substituted by Al and Ga cations [2,3] . With the rapid advances of 
technology, human beings begin to develop devices working in the terahertz (THz) spectrum (0.1 
THz - 10 THz) of electromagnetic wave. It is well accepted that most of naturally synthesized 
magnetic materials can not play roles in the THz domain due to the limitation of their much lower 
resonance frequency. This is the reason that the metamaterials are frequently designed to obtain 
THz permeability [4-6]. Is it really impossible to employ the naturally synthesized materials for 
THz applications? In this paper, we will demonstrate how a very strong hard ferromagnetic L1o- 
FePt (hereinafter called FePt) nanowire can exhibit obvious THz permeability spectra. FePt with a 
tetragonal structure is one of the promising materials with an ultra large areal storage density (^ 
400 Gb/cm?) in the next generation perpendicular magnetic recording industry [7]. Llo-FePt has a 
very large magnetocrystalline anisotropy constants (K) and a very small critical diameter of 
forming single domain without superparamagnetism effect [8-10]. In order to increase the natural 
resonance frequency up to THz, strong shape anisotropy field could add up with the 
magnetocrystalline anisotropy field (H;). Nanowire is a typical shape having strong shape 
anisotropy field. The following table gives the K values and coercivity values (H.) of several 
typical hard magnet and the ease to fabricate them in shape of nanowire by electrochemical 


deposition techniques. 


Table 1: magnetocrystalline anisotropy (K), coercivity (He) , resonance frequency (f;) and ease in 
shape of nanowires !!!! 


Materials K(J/m?) (5! H, (kOe) Ms (A/m) Ease of NWs f. (GHz) Refs. 


Nd:Fe;4B 4.3x10° 73 1.27 x10° No -- al 
SmCo; 1.7x107 400 8.36 x105 No 1132" ist 
BaFe;?0;o 3.2x10° 17 3.82 x10° No 36 Isl 
L1,-FePt 6.6x10° 116 1.14 x10* Good 348" 15] 
MnBi 1.6x10° -- 7.16 x10° No -- [12] 
Ll1,FeCo 1.0x107 -- 1.91x10° Good -- 113] 
L1,-FeNi 1.3x10° -- 1.17 x10° No -- 114] 
(meteorite) 


** patented and to be published. 
2. Micromagnetics theory and simulations 

Micromagnetics approach uses a continuum approximation to analyze magnetization 
behaviors (hysteresis loops), domain structures or dynamic responses of magnetization. 
Micromagnetics describe magnetization on a significant length scale which is large enough to 
replace atomic magnetic moments by a continuous function of position and small enough to reveal 
the transitions between magnetic domains. The micromagnetics simulations are performed using a 
three-dimensional object-oriented micromagnetics framework (OOMMF) by solving the Landau- 


Lifshitz-Gilbert (LLG) equation as a function of time [15]. 


SM, -yM,xH, M,x(M,xH,) 
dt M (1) 


s 


Hey =Hr + Hsn (2) 


, where M, is the magnetization and Her is the effective field taking into account the exchange, 
demagnetization field, anisotropic field and applied field terms. yis the Gilbert gyromagnetic ratio 
(2.21x10° m(A.S) ) and g is the damping constant. In this paper, œ is set as 0.02 for the 
permeability simulations, and is set as 0.5 for the hysteresis loops simulations to save calculation 


time. The default parameters in OOMME for an isolated FePt nanowire are taken in our 


simulations: saturation magnetization Ms is 1.14 x10* A/m, exchange stiffness constant A is 
1x10" J/m, magnetocrystalline anisotropy constant Kı is 6.6x10° J/m? [15] The 
magnetocrystalline anisotropy easy axis is along [001] direction which is set parallel to the wire 
length. The initial magnetization is magnetized along the length direction of nanowire. The shape 


anisotropy constant Ks can be found as: 


K, = 5 HoANM? (3) 

„where AN is 0.5 for a long cylindrical nanowire and is the demagnetization factor difference 
between the hard axis and the easy axis. The anisotropy field due to K; or Ks can be calculated as: 
Hı = 2K/(uMs). A cubic cell size of 2nm x 2nm x 2nm is adopted for all simulations. For 
comparisons, an isolated Fe nanowire is also simulated. The physical parameters for Fe are listed 
as: Ms = 1.7 x10° A/m, A =2.1 x10"! J/m, Kı = 4.8x10* J/m?. Both FePt and Fe nanowires are 100 
nm long with a diameter of 10 nm to ensure strong shape anisotropy. To obtain the permeability 
spectra, we adopt the following procedures. Firstly, a remanent state is chosen as the equilibrium 
configuration of magnetizations, which is obtained by reducing the applied saturation magnetic 
field to zero. Secondly, a weak pulse magnetic field (h(t) = 50 exp(-10"t), where £ in seconds, h in 
A/m) is applied as a small perturbation to the equilibrium state. A(t) is perpendicular to the length 
direction of NWs. The dynamic responses of magnetization in time domain are recorded under 
perturbations. Finally, the fast Fourier transform (FFT) approach is used to convert the recorded 
data into the data in frequency domain, and then the permeability spectra can be drawn. For the 
hysteresis loops simulations, the magnetic field (H) is applied along the length and vary between 
-200 kOe and +200 kOe for FePt nanowire, and vary between -40 kOe and + 40 kOe for Fe 
nanowire. The step size of H field is 0.1 kOe. 


3. Results and discussions 
Firstly, we have simulated the hysteresis loops, as shown in Fig. 1, of single FePt nanowire and 


single Fe nanowire respectively. For these simulations, the external magnetic field is applied along 
the longitudinal direction. Clearly, both nanowires exhibit perfectly rectangular loops and indicate 
both have a strong shape anisotropy field. The easy magnetization direction is along the length of 


nanowire. The remanence (Mr/Ms) is equal to 1.0 (see “A” point), which will be a good 


E) 


equilibrium state for dynamic simulation (i.e. p ~ f ) because almost all magnetic moments 
directions are along the same direction which will give rise to a strong response. The simulated 
coercivity value is 122.51 kOe for FePt nanowire which is larger than the reported value (23 kOe) 
found in the bulk materials with multi magnetic domains and with very weak shape anisotropy 
field [16]. Same fact found in Fe nanowire: the simulated Hc is 10.65 kOe, which is much larger 
than the measured value (356 Oe) of Fe nanowire array with millions of nanowires. Detailed 


explanations were given in our previous paper [17,18]. 
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Fig.1 Hysteresis loops of single FePt and Fe nanowire 


Secondly, the remanent state of FePt or Fe nanowire is illustrated in Fig. 2. Clearly, each 
magnetization vector is pointing toward the longitudinal direction. No domain wall is found. For 
magnetic materials with multi domains, the frequency dependence of permeability will be 
dominated by the domain wall resonance frequency which is far below the natural resonance 
frequency. For instance, we have studied the resonance frequency of domain wall (fow) and natural 
resonance frequency (f.) for NiZn spinel ferrite, and have revealed that fpwis lower than f. by an 
order of magnitude [19].Therefore, it can be inferred that single magnetic domain structure is 
prerequisite to obtain obvious THz permeability for naturally synthesized magnetic materials with 


extremely large K values shown in Table 1. 


Fig. 2 Magnetization equilibrium configurations at remanent state of FePt NW 


When the remanent state (an equilibrium state) show in Fig.2 is used to study the dynamic 


magnetic responses, the complex susceptibility (y=y’- jy”) can be defined as: 


mM, 


h, 0) 


XX) = 


Accordingly, the respective expression for the real parts (u’= 1+y’) and imaginary parts (u” = y^ 
of permeability are given as below [20]. 
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When the following condition is satisfied,the imaginary part of permeability (j^) reaches a 

maximal value. When it happens, this physical phenomenon is termed as “natural resonance". 
Q-—0,-2zf,-yH, — (8) 


The natural resonance frequency (fr) can be calculated according to Eq.(8), in which HO is 


expressed as follows [20]. 


E a 
H,- ^") =H ,(li+—— 
d y a GM? 


1/2 
) x Ha 


—0) 


Since the damping constant (q) is extremely small compared to the product of y and M, see above 
Eqs. 9. H0 can be approximately represented by Heff. In other words, the natural resonance 
frequency (fr) and the permeability spectrum can NOT be adjusted by controlling the damping 
constant, and only be possibly changed by controlling the effective anisotropic field ( Heff) which 
can arise from Hk, Hshape, and (or) other sources, such as stress induced anisotropic field. The 
frequency dependence spectra of permeability are compared for both nanowires and shown in 


Fig.3a. Obviously, the natural resonance frequency of Fe nanowire is about 31 GHz and much 


smaller than that of FePt nanowire (348 GHz = 0.348 THz). According to Eqs.(2,8 and 9), the 
calculated fr value is 31.49 GHz and 0.344 THz for Fe and FePt nanowires respectively, which 
agree very well with the simulated values and corroborate that these observed peaks originating 
from the natural resonance phenomenon related to their Heff field. However, the shape anisotropy 
field (Hsh) plays an important role in enhancing fr value of Fe nanowire, but it can be neglected 
in FePt nanowire because Hsh is much smaller than the magnetocrystalline anisotropy field (Hk). 
Therefore, only one natural resonance frequency is found in FePt. The other resonance peak found 
at 18 GHz is thought as the “edge mode". As indicated in Eq.(7), D is always a positive value. 
Hence, j"is always a positive value (see Eq.6) since damping constant (æ) is generally thought as 
a positive value to represent the energy dissipation during precession. Such a widely accepted fact 
(u" 70) is found in the permeability spectrum of Fe nanowire as shown by the inset in Fig. 3a. 
However, the THz permeability of FePt shows something different: y” firstly reaches a negative 
dip (u" <0) and then jumps back to a normal positive peak at resonance frequency (fr). To 
interpret this abnormal behavior, we need to compare the precessions between Fe and FePt 
nanowires. As shown in Fig.3b, time dependence of Mz component of Fe nanowire 
“symmetrically” oscillates around the zero value, and finally reaches to a static state (Mz=0) when 
the precession of Ms completely stops. For FePt nanowire, Mz firstly shows an “asymmetrical” 
oscillation above zero value and then returns back to a normal symmetrical oscillation around 
zero. Putting the time dependence of Mz and My together will be more useful to demonstrate the 
abnormal precession of FePt nanowire, as shown in Fig.4a. For Fe nanowire, both My(t) and Mz(t) 
follow the same precession behavior as shown by their overlapped time dependence curves in 


Fig.4b. For FePt nanowire, My(t) and Mz(t) are not overlapped at the beginning of precession. 
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Fig.3 THz permeability of FePt nanowire and damping of magnetization component compared 


with Fe nanowire. 
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Fig.4 Damping behaviors of magnetization components of FePt and Fe nanowires. 


By comparing the precession patterns, the THz permeability of FePt nanowire is anticipated to be 


different from those of Fe nanowire. Under the excitation of external pulse magnetic field, the 


magnetization (M) will deflect from its equilibrium direction with an angle (termed as “precession 


angle") and start to precess around this direction (called Lamor precession). In our case, the 


effective anisotropy field direction is along the “x” axis direction. The other one torque is termed 
as “damping torque”, which is pointing toward to the circumference center of precession, and 
reduce the precession angle when the absorbed energy from external field is dissipated away. In 
other word, gradually decreasing precession angle denotes the positive damping constant, vice 
versa. The imagninary parts (j^) of permeability representing the energy loss is therefore always 
positive for normal precession. From Eq. 6, the positive damping constant ensures the positive j^ 


values. Therefore, it is widely accepted as a natural law that jj" can NOT be negative for general 


physical scenarios. However, if the precession angle can be enlarged during an abnormal THz 
precession as shown in FePt nanowire, it is equivalent to have a negative damping constant and 
therefore it will result in a negative u”. Negative jj" phenomenon has also been proposed by us in 
other paper for normal precession where the spin transfer torque of spin-poloarized current were 
used to enlarge the precession angle [21]. Furthermore, from the standing point of energy 
dissipation, when the external energy (e.g. spin transfer torque effect) 1s larger than the dissipated 
energy, it will result in an enlarged precession angle and an equivalently negative damping 
constant. Accordingly, negative jj" will be observed as per Eq.6. In addition, it should be noticed 
that the damping time of magnetization shown in Fig.4 is different for FePt and Fe nanowires. It is 


almost 10 times longer for Fe nanowire. 
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Fig. 5 Variation of total Gibbs energy under perturbation of external magnetic field 
The differences in precession between FePt and Fe nanowire also can be viewed from the 


perspective of Gibbs energy (Ey), as shown in Fig. 5. When the normal nature resonance occurs, 
Ew has a sudden jump due to the energy coupled from the pulse magnetic field excitation, then 
Ew; is dissipated by magnetic loss and the precession angle gradually decreases to zero. Hence, it 
gives rise to positive yy” values. However, several peaks of Ey; are found for FePt, which means 
the precession angle is not gradually decreasing. At some moments, the angle is enlarged which 
can be inferred from the suddenly increased Ew. As stated before, it is equivalent to have a 
negative damping (œ) constant which will result in negative yz” values. Obviously, from the point 


“D” to the point “F” in Fig.5a, Ew: gradually dissipates as normal, therefore, we observe a normal 


positive u” values again. 


4. Conclusions 
Under the remanent state (Mr/Ms) with 1.0, L1lo-FePt nanowire with a large H, value and a 


square M-H loop can exhibit an obvious THz permeability spectrum (p ~ f) with a natural 
resonance frequency (f;) of 0.348 THz. Below some frequency around f., u” is surprisingly found 
negative which is thought due to the abnormal precession behaviors and has been compared with 
an isolated Fe nanowire which shows normal precession behaviors. We believe that if the 
precession angle is enlarged during the Lamor precession, and it is equivalent to have a negative 


damping constant and therefore it will result in negative y” values. 
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